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ABSTRACT 


A  fiber  optics  particle -sizing  system  is  discussed  with  respect  to 
theory  of  operation  and  data  acquisition  and  reduction  techniques.  The 
system  utilizes  a  shadowgraphic  technique  to  determine  the  dimensions 
and  numbers  of  particles  moving  in  a  flow  field.  The  system  is  digital 
in  nature.  Particles  pass  through  a  collimated  laser  beam  and  are 
imaged  onto  a  linear  array  by  a  coaxial  lens.  The  array  is  composed 
of  the  exposed  ends  of  a  fiber  optics  bundle  which  serves  as  a  "link" 
between  the  array  plane  and  the  sensing  photo-detector  modules. 

Being  an  imaging  device,  it  can  measure  a  wide  range  of  particle  sizes 
through  the  proper  selection  of  optics.  Sizes  ranging  from  2  to  1500  fun 
have  been  measured  in  the  course  of  this  project.  Comparisons  of 
holographic  data  of  a  liquid  rocket  injector  and  of  water  spray  nozzles 
are  made  with  the  fiber  optics  system. 
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SECTION  I 
INTRODUCTION 


1.1  PROGRAM  OBJECTIVE 

The  objective  of  this  program  is  to  provide  system  evaluation  of  a 
particulate-monitoring  system  selected  under  Project  BB513  (Particle 
Contamination  Sensor  in  High-Temperature  Regeneratively  Heated  Wind 
Tunnels).  The  program  objectives  called  for  the  monitoring  of  the  size 
distribution  and  particle  flux  of  ceramic  contaminants  emitted  from  a 
large,  ceramic-filled  storage  heater,  with  specific  application  to  the 
Aerodynamic  Propulsion  Test  Unit  (APTU).  These  contaminants  are 
produced  from  the  heat  exchanger  materials  by  thermal  shock  and  by 
the  mechanical  abrasion  which  occurs  during  the  heating  and  cooling 
cycles.  Particulate  contaminants  in  the  test  flow  could  cause  surface 
damage  to  a  test  model  and  also  be  a  source  of  test  data  error.  A  con¬ 
tinuous  contamination  monitor  could  provide  contamination  level  infor¬ 
mation  necessary  for  correcting  the  test  data.  Such  a  system  could  also 
provide  heater  bed  information  and  perhaps  aid  in  the  forecast  of  im¬ 
pending  structural  failure  of  the  bed  by  noting  changes  in  the  dusting 
rate  or  size  distribution.  If  detected  in  time,  some  deterioration  can 
be  repaired  by  heat  treatment. 

The  initial  program  was  subdivided  into  the  five  following  tasks: 

1.  Determine  the  approximate  size  and  quantity  of  contami¬ 
nants  expected  from  such  heaters. 

2.  Determine  the  particle  size  and  concentration  that  would 
be  detrimental  to  wind  tunnel  operations. 

3.  Describe  the  environment  in  which  the  device  will 
operate. 

4.  Survey  the  available  devices,  or  techniques,  for  this 
purpose  and  choose  the  one  most  suitable. 

5.  Evaluate  the  system  operation. 

Under  Project  BF5135  the  first  four  tasks  were  completed  and  the 
fiber  optics  system  (FOS)  was  finally  chosen  as  providing  immediate 
promise  for  application  under  the  unfavorable  operation  conditions  to  be 
encountered.  The  system  had  to  operate  in  an  extreme  temperature, 
velocity,  and  pressure  environment  which  was  made  more  untenable  by 
the  presence  of  density  gradients. 
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12  ELECTRO-OPTICAL  TECHNIQUE 

The  electro-optical  technique  was  originally  developed  and  reported 
(Ref.  1)  by  personnel  at  the  National  Center  for  Atmospheric  Research, 
Boulder,  Colorado.  The  device  was  used  for  measurements  of  size  dis¬ 
tributions  of  cloud  and  precipitation  particles  in  the  range  from 
5  to  1000  Mm  at  low,  subsonic  aircraft  speeds.  The  speeds  expected  in 
tunnel  application  were  an  order  of  magnitude  higher,  while  the  same 
size  range  was  to  be  measured.  For  this  reason  Transistor-Transitor 
Logic  (TTL)  electronics  and  photomultiplier  (PM)  tubes  were  used  for 
the  detector  system.  At  the  time,  this  combination  promised  the  fastest 
response,  which  was  approximately  30  MHz.  A  detailed  description  of 
the  system  electronics  operation  is  outlined  in  Section  5. 1,  Theory  of 
Operation. 

This  report  will  be  concerned  with  the  application  of  the  basic 
instrument  to  tests  at  the  Arnold  Engineering  Development  Center 
(AEDC).  Emphasis  was  placed  on  system  performance  in  typical  wind 
tunnel  conditions.  Data  reduction  methods  were  developed,  correcting 
for  such  factors  as  depth  of  field,  effective  array  width,  average 
particle  velocity,  and  counting  time.  Finally  the  results  of  two  typical 
applications,  a  water  spray  droplet  study  and  a  Propulsion  Wind  Tunnel 
Facility  (PWT)  Dust  Erosion  Study,  including  instrumentation,  will  be 
discussed. 


SECTION  II 
PARTICLE  SAMPLING 


In  general,  particle -sampling  techniques  perform  spatial  and  tem¬ 
poral  averages  of  the  particulates.  All  techniques  will  sample  over  a 
region  in  space  and  time.  Transient  events,  if  they  are  to  be  time 
resolved,  require  fast  sampling  rates.  A  large  number  of  repetitions 
of  the  event  or  large  sample  volumes  are  required  to  obtain  significant 
counting  statistics.  If  the  process  is  quasi  continuous,  sufficiently 
good  time  averages  may  be  obtained  over  a  relatively  small  sample 
area.  APTU  and  other  continuous  tunnels  fall  into  this  category.  The 
FOS  is  eminently  suited  for  this  type  of  application,  as  it  provides 
rapid  counting  over  an  area  and  provides  a  continuous  update  of  particle 
size  information. 
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For  an  understanding  of  the  problem,  the  following  two  basic  cases 
shall  be  discussed:  (I)  a  spatially  averaged  particle  size  distribution 
obtained  in  one  instant  of  time  and  (II)  a  time  average  distribution 
taken  as  average  over  an  area  in  space.  The  difference  between  the  two 
distributions  is  considerable  under  certain  conditions. 

Case  I  deals  with  scanning  by  photography,  holography,  and  other 
short-time  procedures  (Ref.  2).  These,  techniques  record  a  volume  of 
the  flow  in  a  period  of  time  which  is  short  with  respect  to  the  character¬ 
istic  time  of  the  particles,  t,  to  prevent  blurring  (Fig.  1,  Appendix  I). 

D 


where  D  =  particle  diameter  and  V  =  particle  velocity.  The  distribution 
obtained  from  such  a  procedure  is  density  weighted  (i.  e. ,  the  distribution 
function  gives  just  the  average  number  density  of  particles  of  a  particle 
size). 


n(D) 


number 

volume 


Case  II  requires  real-time  particle  analyzers,  including  the  FOS. 
Although  for  analysis  purposes  the  latter  contains  also  a  sampling 
volume,  one  can  assume  that  the  volume  is  replaced  by  an  area  through 
which  the  particle  flux  streams  (Fig.  2).  Each  partide_falls  into  a 
velocity  range  Vj  =  ±  AVj  and  a  diameter  range  Dj  =  ±  ADj 

associated  with  it.  As  each  particle  is  swept  into  the  active  area,  it  is 
sized  by  some  means.  The  resulting  distribution  is  weighted  by  the 
particle  flux  as  follows: 

number 

j.(D)  =  - 

area-sec 


The  two  distributions  are  not  independent,  as  they  are  related  by 
the  continuity  equation.  If  the  volumes  are  assumed  to  be  infinitesi¬ 
mally  small,  the  continuity  equation  is  written  as 


div  nv  =  sr  ^  4 

at  tPy 


24 


Note  that  n  and  V  are  functions  of  D.  Assuming  steady-state  particle 
fields,  the  continuity  equation  reduces  to 
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or 


nV  =  J  =  Constant 


Therefore,  the  number  density  of  case  I  systems  is  related  to  the 
particle  flux  of  case  II  systems  by  the  particle  velocity.  In  most  appli¬ 
cations  where  the  flow  is  obtained  by  the  acceleration  of  a  gas  through 
a  nozzle  or  around  a  model,  one  can  not  expect  to  have  an  equilibrium 
between  particle  and  gas  velocity  fields.  Since  drag  is  a  surface  effect 
and  inertial  forces  are  volume  dependent,  a  velocity  distribution  would 
be  expected  for  the  variously  sized  particles.  In  general,  a  relative 
particle  size  distribution  obtained  by  a  case  I  system  will  be  different 
from  one  obtained  by  a  case  II  system,  for  accelerated  particle  systems. 


2.1  GENERALIZED  COUNTING  -  A  STATISTICAL  APPROACH 

The  generalized  quantity  which  is  desired  to  characterize  the  pro¬ 
cess  is  the  distribution  function  based  on  particle  flux.  The  total  mass 
flow  can  be  obtained  from  the  distribution.  This  quantity  is  the  number 
of  particles  passing  through  a  unit  area  in  a  unit  time  for  a  unit  range 
of  particle  diameters.  The  FOS  automatically  detects  these  particles, 
sizes  them,  and  adds  them  to  the  appropriate  channel.  A  data  reduction 
scheme  will  be  described  which  was  used  to  determine  the  normalizing 
factors  for  each  channel.  That  is,  the  number  stored  in  the  ith  channel 
is  representative  of  a  particular  sensing  area  and  size.  The  exact  area 
must  be  calculated  to  form  the  desired  distribution  function.  To  do  this, 
the  probability  of  a  particle  of  a  given  size  and  at  a  given  location  on  the 
conjugate  array  plane  being  detected  and  sorted  into  a  given  channel  will 
be  calculated.  Once  this  function  has  been  determined,  average  quanti¬ 
ties  (e.  g. ,  the  average  particle  size  within  a  channel,  or  the  effective 
array  length)  can  be  calculated. 

In  Fig.  3  particles  are  shown  to  be  moving  orthogonal  to  the  fiber 
array  in  the  y  direction.  They  are  illuminated  from  the  z  direction  by 
collimated  light.  The  motion  will  carry  the  shadows  past  the  array 
located  symmetrically  along  the  x  axis.  Whether  or  not  a  fiber  optic 
module  is  triggered  is  determined  only  by  the  intensity  in  the  shadow 
region  and  by  its  x  position.  If  we  assume  no  diffraction  effects,  then 
the  shadow  will  be  perfectly  black,  and  geometrical  considerations 
will  adequately  describe  the  system  response. 
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2.2  PROBABILITY  MODEL 

The  purpose  of  the  probability  model,  which  is  to  be  developed,  is 
to  provide  a  method  for  determining  the  normalizing  factors  for  each 
channel;  this  method  depends  on  the  effective  array  length,  the  sensing 
volume  depth,  and  the  mean  particle  diameter.  The  generality  of  the 
model  allows  for  more  complex  intensity  patterns  in  the  array  plane  to 
be  analyzed.  Consider  a  particle  of  diameter  D  and  position  x  in  the 
array.  Pj  (D,  x)  is  defined  as  the  probability  of  a  particle  of  diameter  D 
producing  a  count  in  the  i**1  channel  when  located  at  position  x.  P  can 
talce  on  values  only  of  0  or  1  {no  or  yes).  To  produce  a  count  on  the  i**1 
channel  the  particle  presence  must  shadow  and  trigger  exactly  i-fiber 
optic  modules.  If  it  does,  P  =  1;  otherwise,  P  =  0.  The  density  of 
particles  will  be  assumed  to  be  uniform  over  that  segment  of  the 
x  axis  containing  the  array  {i.  e. ,  -X/2  <  x  <  X/2).  P  is  also  a  function 
of  threshold  level  in  a  manner  which  will  be  shown  later. 


2.3  DISTRIBUTION  FUNCTION 


The  desired  distribution  function  is  discrete  in  nature  (i.  e. ,  there 
is  one  value  for  each  active  sizing  channel).  However,  one  may  assume 
that  the  discrete  values  are  points  taken  from  a  continuous  function. 


dx  H  7.  dD  dt 


In  the  case  of  limited  depth  of  field,  where  cutoff  is  provided 
mechanically  either  by  physical  masking  or  by  electro -optical  coinci¬ 
dence  counting  based  on  scattered  radiation,  the  z -dependence  can  be 
integrated  out  to  provide 


d3N 


dx  di)  dt 


1 

*  h 


where  =  effective  depth  of  field. 

If  a  time  average  is  performed  with  the  assumption  that  no  change 
in  particle  flux  is  incurred  over  the  sample  time,  then 


r.  _  J*. 

>  dx  dl) 


1 


where  T  =  sample  time. 
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In  a  similar  manner,  if  Pj(x,  D)  can  be  determined,  the  effective 
diameter  and  effective  array  width  can  be  also  obtained.  This  can  be 
done  by  forming 

I;  =  /  Pi(x,D)dx  dD 


For  a  symmetric  array  of  7  fibers,  y-2  of  which  are  active. 


(y-2)b 

Ad 

2  P,(x£.D„)4»4D 

m=0 


The  distribution  function  is  now  reduced  to  measurable  and  calcu¬ 
lable  quantities  where  Nj  is  the  stored  counts  in  the  ith  channel  after  a 
period  of  time,  T.  The  above  equation  provides  the  basis  for  a  data 
correction  computer  program. 

In  a  similar  manner  Dj  can  be  calculated 

(y-2)b  (y— 2)b 

2Ax  Ad 

2  d  l0  D- 


i, 


(y-2)b 
2Ax 
2  1 
1=0 


f-  = 


I;L.T 


This  is  just  the  average  size  of  a  particle  which  would  contribute  to  the 
counts  in  the  i**1  channel.  In  general  (see  Fig.  4), 

U.  =  ib  =  iMDf 

where  b  =  fiber  size  in  the  object  space,  M  =  system  magnification, 

Df  =  fiber  diameter.  The  effective  array  width  is,  similarly. 


_yb_ 

2  Ax 

w<  '  24xei  p*(xe’r’i) 


In  general  this  quantity  is  approximately  (30-i)b  (see  Fig.  5). 
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To  this  point,  no  mention  has  been  made  of  how  Pi(x,  D)  is  calcu¬ 
lated.  A  technique  will  now  be  developed.  Consider  the  fiber  array 
with  a  particle  located  at  some  point  x  as  measured  from  the  array 
center  (Fig.  6).  To  determine  whether  a  fiber  will  trigger  or  not,  one 
need  only  consider  the  case  where  the  number  of  fibers  shadowed  is  at 
maximum  during  a  particle  passage.  This  will  occur  when  the  particle 
center  lies  on  the  array  axis.  All  fibers  lying  totally  within  the  particle 
shadow  will  trigger  at  some  time  during  the  particle  transit.  These 
fibers  partially  occulted  must  be  examined  in  more  detail. 

Figure  7  illustrates  the  area  of  overlap  for  a  round  particle  and  a 
round  fiber.  If  the  area  of  overlap  is  greater  than  that  minimum  area 
required  for  triggering,  the  fiber  will  count;  otherwise,  it  will  not. 

This  argument  is  strictly  true  only  for  a  perfectly  black  shadow. 

D  ±  2X 

The  edges  of  a  particle  are  located  at  d±  =  — ^ —  fiber  diameters. 

The  plus  sign  corresponds  to  the  edge  in  the  positive  x  direction  from 
the  particle  center,  whereas  the  negative  sign  corresponds  to  the  edge 
in  the  negative  x  direction. 


The  term  d±  can  be  expressed  in  terms  of  an  integer,  J^.,  and  a 
positive  quantity  less  than  unity.  d±  can  be  expressed  as  Jf+e-t,  0  <e.fc<  1. 
The  integer  terms  give  the  total  number  of  fibers  completely  shadowed 
for  an  indeterminately  long  array,  whereas  the  area  of  overlap  must 
be  calculated  based  on  to  determine  if  an  additional  one  or  two  fibers 
will  contribute.  In  practice,  since  the  array  is  limited  in  extent,  the 
number  of  fibers  occulted  must  be  checked  to  see  if  it  exceeds  the 
active  array  length.  If  it  does,  no  counts  will  be  registered  since  the 
end  fiber  reject  mechanism  will  not  allow  the  particle  to  be  sized. 


From  geometrical  considerations  only,  the  total  area  of  overlap 
(Fig.  7)  is  given  by 


b20, 


A  = 


Aj  +  A  2 


0,  =  cos"I(2A]/b) 

0, 


D20, 

-  AAi  +  v 


—  AA2 


cos-](2A2/D) 


-A„  = 


D2  -  b2 


l/2[f  -  l/2(b  +  D)]  + 

-A]  =  l/2[f-l/2(b  +  D)]  +  5Lc_“1/-b".  m 


8le-  l/2(b  +  D)], 
b2  -  D? 
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where 


A 


If  e±  and  d-j.  are  substituted  into  the  equation  for  A,  and  if  b  =  1, 
the  result  is  area  A±. 

A  new  variable,  K-j.,  will  be  now  introduced,  which  is  defined  as 
follows: 


A+ 

K+  =  0  if -  <  1  -  thresh 

1  *<1)2  ~ 


+ 

K+  =  1  if  -  >  ]  -  thresh 


where  "thresh"  is  the  fraction  of  the  average  intensity  which  the  instan¬ 
taneous  intensity  must  fall  in  order  to  trigger  the  electronics.  This 
merely  states  that  will  take  on  a  value  of  unity  if  threshold  triggering 
conditions  are  met.  Otherwise,  it  is  zero.  Now  P^  can  be  determined. 

Pj(x,D)  =  0  if  J+  +  K+  or  J_  +  K_  >  15  or  (J+  +  J_  +  J_  +  K_)  #  i 

P;(x,D)  =  1  if  J+  +  K+  and  J_  +  K_  <  15  and  (J+  +  K_  +  J_  +  K_)  =  i 

can  now  be  substituted  back  into  the  equations  for  Dj  and  W^.  A  com¬ 
puter  program  for  performing  this  calculation  and  correcting  the  data 
and  a  sample  computer  printout  with  accompanying  explanation  are 
shown  in  Tables  I  and  II,  respectively  (Appendix  II),  in  addition  to  the 
average  size  in  the  i  channel  and  the  effective  array  width. 

The  average  bin  size  (i.  e. ,  that  range  of  size  over  which  particles 
will  be  sized)  is  in  the  i ^  channel  and  is  given  by  the  difference  between 
succeeding  mean  channel  sizes. 


b 


m+Ej-Dj),  i  =  i 

1  <  i  <  30 

K(b  +  —  ^9),  i  —  30 
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Calculations  of  effective  sample  area,  distribution  function,  parti¬ 
cle  flux,  and  particle  density  must  be  deferred  until  a  complete  discus¬ 
sion  of  the  depth  of  field  has  been  made. 


2.4  LINEAR  APPROXIMATION  TO  DETECTION  PROBABILITY 

An  approximate  solution  will  now  be  presented.  In  this  case  the 
assumption  is  that  the  array  is  truly  linear  in  the  sense  that  the  fibers 
are  b  units  wide  in  the  x  direction  and  are  infinitesimal  in  the  y  direc¬ 
tion.  Therefore,  an  overlap  of  a  particle  will  produce  a  reduction  in 
intensity  proportional  to  the  fraction  of  the  fiber  covered  in  the  x  di¬ 
rection  only  (see  Fig.  8). 

If  i  fibers  have  been  triggered,  the  range  of  particle  sizes  which 
could  have  been  responsible  may  be  determined  by 

D-,  max  =  b[i  +  2(1  -  thresh)] 

Di  min  =  b['  >  (2  thres^ 


The  maximum  occurs  when  the  particle  totally  covers  i  fibers  and 
covers  just  less  than  the  amount  required  to  trigger  each  of  the  two  ad¬ 
jacent  fibers.  In  a  similar  manner,  the  minimum  detectable  size 
occurs  when  a  particle  covers  i  fibers  but  covers  the  outermost  two 
fibers  just  enough  to  cause  them  to  trigger. 

The  size  range  of  particles  as  sized  in  the  itlri  channel  is 
13;  -  1  <  D.  <  Dj  ^  1  where  Dt  =  — maX I  ^  min 

Of  course,  the  probability  of  a  particle  at  the  extreme  values  being 
counted  declines  to  a  vanishingly  small  value.  The  actual  probability 
distribution  can  be  obtained  from  the  previous  analysis.  The  total  dis¬ 
tance  over  the  array  which  a  particle  of  size  D  can  move  and  trigger 
is  given  by 


£,(D)  =  /~p.(x,D)dx 

— oc 


Figure  9  shows  a  plot  of  the  relative  probability  of  triggering  versus 
particle  diameter  for  the  linear  approximation. 
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SECTION  III 

ELECTRO-OPTICAL  REQUIREMENTS 


3.1  APERTURE  REQUIREMENTS 


3.1.1  Minimum  Numerical  Aperture 


The  resolution  of  an  imaging  system  in  terms  of  resolvable  spa¬ 
tial  frequency  is  proportional  to  the  available  numerical  aperture. 
Therefore,  in  system  design  attention  must  be  given  to  the  size  and 
location  of  the  lens  so  that  adequate  collection  angles  are  provided. 


A  system  must  have  a  minimum  numerical  aperture  (N.  A. ) 
which  can  be  calculated  based  on  the  object  geometry  (Ref.  3).  For 
a  grating-like  object  the  limit  of  observation  is 


0.61 


A 

n  sin  d 


0.61  A 
N.  A. 


For  incoherent  radiation  and  for  coherent  radiation,  it  is 

.  ,  0.77A  _  0.77A 

''  n  sin  6  N.  A. 

In  both  cases  a  circular  aperture  is  assumed.  While  the  constant  fac¬ 
tors  are  somewhat  arbitrary,  as  they  depend  on  the  form  of  the  object, 
the  aperture,  and  the  sensitivity  of  the  receptor,  they  are  commonly 
used  as  the  resolution  limit. 


3.1.2  Depth  of  Field 

A  particle  illuminated  by  uniform  plane  radiation  produces  a  dif¬ 
fraction  pattern  at  a  distance  z  away  characterized  by  the  particle  size 
and  the  wavelength.  The  scaling  parameter  is  the  farfield  and  is  given 
by  D2  /  A. . 

At  the  distance  of  approximately  one  farfield,  the  intensity  pattern 
still  resembles  the  geometric  shadow.  Beyond  this  point  diffraction 
effects  are  dominant.  Therefore,  the  depth  of  field  for  a  particle  would 
be  expected  to  be  some  constant  multiple  of  this  characteristic  distance. 
Figure  10  shows  a  particle  of  diameter  D  imaged  with  a  lateral  magnifi¬ 
cation  M.  For  small  distances  about  the  image  point  the  longitudinal 
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magnification  is  M  .  Assuming  a  perfect  imaging  system,  the  disturb¬ 
ance  at  one  farfield  from  the  particle  is  transferred  by  the  lens  to  a 
point  approximately  that  distance  away  from  the  particle  image. 

But 


which  is  a  farfield  based  on  the  image  particle  diameter.  The  intro¬ 
duction  of  a  diffraction-limited  lens  system  does  not  effect  changes  in 
the  depth  of  field  as  sensed  at  the  array,  as  the  overall  intensity  in 
the  shadow  region  decreases  with  z.  The  points  of  electronic  cutoff 
will  be  approximately  given  by  z  =  CD2/4A  where  C  is  a  function  of  the 
threshold.  At  the  limits  given  by  the  above  equation  the  intensity,  on 
the  average,  will  have  risen  above  that  value  required  to  trigger  a 
module. 

By  combining  the  numerical  aperture  requirement  with  the  depth 
of  field  relationship  one  obtains  an  expression  which  is  a  measure  of  the 
adequacy  of  the  lens  system.  Refer  again  to  Fig.  10  for  a  simple  lens 
system  with  effective  lens  diameter  JL .  A  quantity  R  will  be  formed 
which  is  the  ratio  of  the  required  numerical  aperture  for  cutoff  to  the 
available  numerical  aperture.  That  is,  for  a  given  position,  z,  a 
particle  of  some  size  D(z)  will  undergo  cutoff  as  is  given  by  the  depth- 
of -field  relation.  This  particle  has  associated  with  it  a  minimum 
numerical  aperture  which  is  given  by  the  expressions  for  the  resolution 
limit,  as  follows: 


Figure  11  shows  a  plot  of  R  as  a  function  of  z  for  A  =  0.  6328  pm 
and  zQ  =  50  cm  and  a  lens  diameter  l  =  5  cm.  The  lens  in  the  example 
provides  adequate  numerical  aperture  for  particles  with  a  cutoff  dis¬ 
tance  of  0.  235  cm.  This  corresponds  to  a  minimum  resolvable  particle 
of  48  Pm.  A  practical  limit  occurs  for  f /  =  1  which  would  correspond 
to  a  particle  diameter  of  1  Pm  in  the  visible  spectra. 

3.2  THE  OUT-OF-FOCUS  OPAQUE  PARTICLE 
3.2.1  Intensity  Distribution  Calculation 

A  method  for  calculating  the  expected  intensity  distribution  on  the 
fiber  optic  array  would  allow  parametric  variation  for  system  design 
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and  analysis.  Important  parameters  are  lens-to-object  distance,  mag¬ 
nification,  source  wavelength,  lens  diameter,  and,  of  course,  particle 
diameter.  For  example,  it  would  be  of  considerable  interest  to  know 
in  advance  if  a  10 -pm  particle  could  be  measured  and  sized  accurately 
in  a  tunnel  where  only  a  5-cm  lens  would  be  available  at  a  distance  of 
20  cm.  It  would  also  be  desirable  to  know  the  distance  over  which  it 
could  be  counted  for  an  accurate  depth-of -field  correction. 

A  brief  development  of  the  equations  will  now  be  given.  Figure  12 
illustrates  the  basic  geometry.  This  consists  of  a  simple  lens  system 
imaging  points  from  a  distance  s  away  from  the  lens  onto  the  array 
plane  located  at  a  distance  z  =Ms,  where  M  is  the  system  magnifica¬ 
tion. 


The  particle  is  assumed  to  be  an  opaque  disk  located  on  the  optical 
axis.  The  disk  of  diameter  D  is  located  a  distance  z  away  from  the  lens 
in  the  object  field  and  can  be  specified  by 

z  =  s  +  A  / 

For  a  typical  application,  Az  would  be  the  distance  from  duct  centerline 
of  a  uniform  particle  flow  field  monitoring  system.  Coordinates  of  the 
disk  plane  are  (rQ,  ©q,  z),  those  of  the  lens  are  (r^  ©j,  0),  and  those 
of  the  array  plane  are  (r2,  ©2,  z').  Coherent  light  of  unit  amplitude 
intensity  is  assumed  to  be  incident  from  the  left. 

Babinet's  principle  is  used  to  determine  the  distribution  of  the  light 
amplitude  at  the  lens.  That  is,  the  light  amplitude  distribution  resulting 
from  a  disk  can  be  computed  from  the  unperturbed  distribution  plus  that 
of  a  hole  in  the  following  manner: 

lfn  =  '■n  +  Lh 

where 

li0  =  wave  amplitude  without  object 
II jE  =  amplitude  due  to  presence  of  hole 

U„  =  amplitude  due  to  presence  of  disk 

On  the  assumption  that  the  lens  is  at  least  a  farfield  from  the  disk, 
the  Fraunhofer  equations  can  be  used  to  propagate  the  disturbance 
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caused  by  the  disk  to  the  lens  (Ref.  4).  The  disturbance  then  becomes 


L'i){r].0j)  =  exp  j  k  z 


exp  (jkrj^/27.) 

jAz 


/D\2  J  ;(i7Drj/Azi 
W  (Dr  j/  2Az) 


where  k  =  2tt /  A.  is  the  propagation  constant. 


The  disturbance  is  propagated  through  a  simple  lens  of  focal  length 
f  and  diameter  £  and  results  in  a  disturbance  amplitude  U^. 


|-n  =  CXP 


UD 


where 


circ 


jl.fi  <  f/2 
lO.  rj  >  P/2 


Since  the  lens  is  less  than  a  farfield  from  the  array  plane,  the  lens 
is  the  characteristic  size,  and  the  Fresnel  equation  must  be  used  to 
propagate  the  disturbance  through  z',  arriving  at 

ro(r,,0o)  =  — J-r  f  f  Un  exp  [r?  +  r2  -  2rjr2  cos  0] 

~  ~  ~  jAr1  »  0 


If  the  modified  Bessel  function  of  the  zeroth  order,  IQ,  is  used, 
this  can  be  simplified  as  follows: 


1  P/2 

U2(r2,02)  =  - -  /  U|'}  exp 

j\  /.  ° 


r  jdr 


1 


where 


2n 


I(1  =  f  exp  [x  sin  0]Aq 

O 

krl  r2 
x  =  -j  — r- 


The  equation  is  solved  by  numerical  integration  and  is  tabulated  in 
Table  HI.  As  the  photomultiplier  tubes  are  square  law  detectors,  the 
important  quantity  is  the  intensity,  which  is  given  by 

I(r2)  =  U2U2  =  (^‘(Re2  +  Im2) 
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where  Re  and  Im  are  the  real  and  imaginary  parts  of  the  amplitude. 
After  some  manipulation  they  can  be  put  in  the  following  form: 

Re  =  /F/2rJJ»(“T-L)  cos  (Ar0dri  -  //2ji(^)Jo(^)  sin(0riK 


lm  -  //2rlJo(-^)sinH)+  T  Jo(^)  cos  (®r])dfj 


where 


■' 4  Kf  - 1) 


The  intensity  distribution  for  a  sample  particle  geometry  is  given 
in  Fig.  13.  This  corresponds  to  the  image  of  a  20-pm  particle  44  cm 
away  from  a  22 -cm  focal  length  lens  of  diameter  5.  5  cm.  The  system 
magnification  is  unity;  therefore,  the  particle  is  in  focus.  The  pres¬ 
ence  of  the  intensity  variations  is  a  result  of  the  finite  lens  diameter. 
The  true  particle  radius  is  indicated,  as  well  as  that  of  50-  and 
30-percent  thresholds.  A  30-percent  threshold  would  undersize  this 
particle  by  approximately  10  percent.  The  strong  central  intensity 
peak  would  cause  certain  of  the  fiber  optic  models  to  indicate  the 
absence  of  a  particle.  However,  as  will  be  explained  later  in  more  - 
detail,  the  gating  system  would  still  size  this  particle  near  the  indi¬ 
cated  points.  Studies  involving  a  television  (TV)  monitoring  system 
illustrate  this  effect. 

3.2.2  Television  Analysis  of  a  Particle  Passage 

For  this  study  a  300 -/um  particle  was  placed  on  a  glass  disk.  It 
was  illuminated  by  a  collimated  helium-neon  (He-Ne)  laser  beam  and 
imaged  onto  a  vidicon  tube.  The  image  was  then  displayed  on  a  TV 
screen,  where  the  image  could  be  photographed  (see  Fig.  14).  One 
of  the  raster  lines  could  be  selected  to  display  the  intensity  profile. 

A  series  of  photographs  and  intensity  profiles  was  taken  in  sequence  to 
simulate  particle  passage.  For  the  sequence  shown,  the  particle  and 
array  were  not  centered  on  the  beam  profile.  While  all  particles  with¬ 
in  the  geometric  shadow  are  not  at  trigger  level  simultaneously,  they 
do  reach  that  value  sometime  during  the  particle  passage. 
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3.3  CUTOFF  LENGTH 

The  points  along  the  optical  axis  which  define  the  limits  of  particle 
detection  are  very  important  in  determining  the  actual  particle  flux  for 
any  particular  size  of  particle.  The  data  reduction  program  includes  a 
semiempirical  equation  which  expresses  the  active  length  as  function  of 
particle  size.  For  the  ideal  system,  where  the  optics  are  not  band 
limited  and  are  aberration-free,  a  reasonable  form  to  assume  is  that 


This  is,  of  course,  a  statement  of  the  scaling  by  the  farfield  relation  as 
previously  shown. 

Attempts  were  made  to  determine  the  actual  depth  of  field  for  each 
particle  size  by  several  means.  This  included  using  a  traverse  system 
to  accurately  move  a  particle  aspirator  through  the  focal  region.  The 
data  were  then  used  to  determine  average  cutoff  lengths  as  a  function  of 
size.  Figure  15  shows  the  traverse  system  used,  including  the  aspirator. 
The  aspirator  consisted  of  an  air  source  consisting  of  3/8-in.  tubing 
coaxial  with  the  1/8 -in.  particle  tube.  Plant  air  at  approximately  150 
psig  was  used.  A  uniform  supply  of  particles  was  provided  by  a  rotating 
hopper/ dispenser  (see  Fig.  16).  A  small  12-v  d-c  motor  rotated  a 
rubber- surfaced  disk  at  approximately  1  rpm.  Particles  were  deposited 
on  the  disk  from  a  tube,  the  rate  being  adjusted  by  the  clearance  between 
the  tube  and  the  surface.  A  pickup  was  provided  diametrically  opposite 
the  dispenser  tube. 

The  experimental  cutoff  was  obtained  by  dispensing  equal  amounts  of 
a  particle  sample  at  various  points  along  the  optical  axis.  The  particle 
distribution  itself  contained  a  broad  range  of  sizes  to  produce  statistically 
significant  particle  counts  in  each  channel.  The  distribution  can  be  used 
to  obtain  the  average  cutoff  length  in  the  following  manner: 

2£  11- (z)  z. 

~n  »  1 


where  n^(z)  is  the  counts  received  in  the  i**1  channel  at  a  position  z.  The 
sum  is  over-all  z.  The  absolute  value  of  the  distance  must  be  used; 
otherwise,  the  first  moment  is  obtained,  which  is  approximately  zero. 
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Figure  17  shows  a  plot  of  cutoff  length  as  a  function  of  particle 
diameter.  A  rather  drastic  departure  from  the  theoretical  dependence 
is  noted  for  the  lower  channels.  Adequate  numerical  aperture  was  pro¬ 
vided,  which  would  rule  out  the  effect  due  to  bandpass  limitations.  The 
fragmentation  of  the  image  caused  by  beam  nonuniformity  may  play  a 
role.  .Furthermore,  the  particles  were  restricted  to  a  region  with  a 
diameter  on  the  order  of  5  mm,  and  no  recirculation  or  fanout  was 
noted  which  could  account  for  the  abnormally  large  depth  of  field. 

That  is,  the  indicated  cutoff  region  is  an  average  value,  as  the  particle 
could  be  located  ±2.  5  mm  about  any  particular  position. 

Electronic  noise  is  a  likely  source  as  noise  events  could  cause 
single  or  multiple  triggers.  A  Poisson  distribution,  which  statistically 
describes  rare  events,  can  be  made  to  fit  the  lower  channel  counts, 
which  is  further  evidence  that  the  counts  are  the  result  of  random 
events.  By  increasing  the  feedback  at  the  comparator,  random  noise 
could  be  made  to  combine  several  single  channel  counts  into  a  higher 
channel  count.  If  the  particles  were  to  remain  in  the  sample  area 
long  enough,  the  count  would  closely  approximate  the  true  particle 
size.  This  would  at  the  same  time  move  the  cutoff  point  further  out. 
Solutions  to  this  problem  will  be  discussed  in  the  conclusion. 


3.4  BEAM  UNIFORMITY 

The  fiber  optic  system  is  by  design  relatively  insensitive  to  illumi¬ 
nation  nonuniformities.  This  is  a  direct  result  of  the  fact  that  each 
fiber  has  its  own  photomultiplier  (PM)  tube  and  discriminator  circuit 
which  can  be  individually  adjusted  to  the  average  intensity  level  at  that 
point.  Also,  changes  in  the  local  light  level  attributable  to  nonuniformity 
of  the  beam  coupled  with  vibration  or  some  other  form  of  beam  drift  are 
self- compensating,  as  the  reference  adjusts  to  slowly  changing  light 
levels.  The  reference  level  is  essentially  provided  by  filtering  of  the 
quiescent  signal.  Therefore,  any  slow  changes  in  level  will  be  auto¬ 
matically  nulled  out  as  the  trigger  level  rises  and  falls  with  the  fluctu¬ 
ation. 

However,  for  the  following  reasons  there  is  a  requirement  for  a 
degree  of  beam  uniformity.  First,  the  self-compensating  characteristics 
just  outlined  have  limits,  especially  when  one  is  dealing  with  low  signal- 
to-noise  ratios  (S/N).  Second  (remembering  that  the  beam  is  nonuniform), 
diffraction  effects  cause  spreading  of  the  light  so  that  diffracted  light  from 
a  point  of  high  intensity  in  the  object  space  may  interfere  with  low -intensity 
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light  in  the  image  space.  This  effect  could  cause  triggering  on  a  dif¬ 
fraction  ring  instead  of  on  the  particle  "edge. 11  Such  a  phenomenon  can 
be  readily  demonstrated.  Figure  18  is  a  photograph  of  a  particle  in  a 
Gaussian  beam.  The  latitude  of  the  film  has  been  exceeded,  giving  the 
impression  of  a  uniform  beam.  The  particle  center  is  not  coincident 
with  the  illuminating  beam,  and  the  particle  is  out  of  focus.  The  pas¬ 
sage  of  this  particle  across  the  array  A  A'  would  probably  leave  many 
of  the  fiber  optics  located  near  the  center  untriggered,  resulting  in  an 
undersized  count. 


If  the  array  were  located  at  BB'  instead,  the  particle  passage 
would  probably  result  in  one  oversize  count  as  indicated  by  the  passage 
of  the  first  diffraction  ring  and  by  a  count  near  that  of  the  actual  parti¬ 
cle  size.  Figures  19  and  20  indicate  intensity  profiles  for  a  particle 
in  uniform  and  Gaussian  illumination,  respectively. 


3.4.1  Uniformity  Criteria  for  Gaussian  Beam 

The  Gaussian  beam  profile  is  easily  obtained  from  most  lasers  and 
possesses  the  unique  property  that  it  retains  its  Gaussian  form  as  it 
propagates. 

At  the  1/e2  point  of  an  initially  collimated  beam  of  diameter  2? 1,  the 
intensity  profile  is  given  by  the  following  equation  (see  Fig.  21): 

,(r-2)  -  Gr)I°exp(-^) 

where 


At  the  focal  plane 

_ 

w°  “  fffl 

where  f  =  focal  length  of  the  lens  and  \  =  wavelength  of  light  (Ref.  5). 

If  a  constraint  is  made  on  the  uniformity  of  the  illumination  for 
variations  in  z  such  that 


I(o,z)  ,  -  i 

— j -  >  a  where  a  <  1 
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then  this  provides  the  following  constraint: 


W 


2 

o 


For  example,  if  z  =  10  cm,  a  =  0.  9,  and  X  =  0.  6328  Mm,  then  the 
f  /  >  642.  If  a  laser  which  has  any  output  diameter  of  1  mm  is  used,  the 
lens  must  have  a  focal  length  greater  than  642  mm.  If  this  were  the 
only  constraint,  the  lens  might  be  avoided  entirely. 

Uniformity  of  illumination  in  the  radial  direction  must  also  be  main¬ 
tained,  as  was  discussed  earlier.  Demanding  that  the  intensity  not  fall 
below  a  fraction,  /3,  of  the  centerline  intensity  results  in  the  following 
constraint: 

2  — 2rm 2  _  -L2 

°  ~  2M2£n/3 


where 


rm  =  radius  of  array  in  object  space 

L  =  length  of  array  in  object  space 
M  =  magnification  of  system 


This  may  also  be  rewritten  in  terms  of  the  effective  f/,  as  before: 

{/  >  1  _  Til.  1 

2A  <-£n/3)K  ”  4  MA  (-fn/3)'1 

In  general,  the  more  demanding  constraint  will  be  that  of  radial 
uniformity,  as  the  two  expressions  are  not  completely  independent. 
Depth  of  field  is  related  to  particle  size  as  the  square  of  the  particle 
diameter.  The  array  width,  2rm,  as  seen  in  the  object  space,  is 
approximately  equal  the  maximum  countable  particle  diameter.  With 
these  relations  the  ratio  of  w0's  can  be  obtained  from  the  a  and  0 
criteria.  With  reasonable  assumptions  concerning  the  values  of  the 
constants  involved,  the  relative  importance  of  these  constants  can  be 
estimated  as  follows: 


wo<0) 


y 

7T  \l-aJ 


(-2Cn/3) 
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Let  C  =  2.  5,  a  =  0.9,  and  j3  =  0.  9;  then 

(v2(a)  =  0.505  w^(^3) 


Because  of  this  relation,  care  need  be  taken  only  in  maintaining 
sufficient  radial  uniformity  in  the  knowledge  that  the  beam  will  be  suffi¬ 
ciently  uniform,  axially.  This  is  especially  true  if  the  cutoff  is  mechani¬ 
cally  limited,  as  the  maximum  depth  of  field  will  never  be  achieved. 

3.4.2  Power  Requirements 


In  system  design,  the  power  requirement  determines  the  laser  to 
be  used.  In  particularly  high-loss  systems,  as  many  wind  tunnel  appli¬ 
cations  can  be,  the  knowledge  of  the  minimum  power  required  provides 
a  basis  for  calculation  of  actual  power  requirements. 

By  integrating  the  intensity  equation  out  to  some  radius  r,  the 
power  contained  within  that  radius  is  obtained  as  follows: 

r  imM2  / 

P(r,o)  =  /  27rrI(r,o)  dr  =  — - —  |l  -  exp 


where  Im  is  the  minimum  intensity  level  which  will  be  encountered  in 
the  array.  The  total  power  is  obtained  by  letting 


lim  P(r,o)  =  P^ 

r-*oc 


l  M2  n w2 
m  o 

~W~  ~ 


Tests  with  various  P's  and  corresponding  system  magnifications  sub¬ 
stantiated  the  above  relationship.  The  equation  does  not  include  the 
effects  of  wavelength  change,  which  is  contained  in  the  expression  for 
Im.  A  knowledge  of  the  phototube's  photo-cathode  characteristics  can 
be  used  to  predict  illuminance  requirements  if  an  experimental  point 
is  obtained.  For  example,  at  0.  6328  Mm,  Im  -2.7  mw/ cm^.  This 
requires  16.  5  mw  of  total  power  for  P  =  0.  7. 

At  this  point  a  comment  should  be  made  about  the  use  of  cylindri¬ 
cal  lenses.  It  is  obvious  that  much  of  the  power  of  a  system  designed 
around  spherical  lenses  is  wasted  since  the  array  is  linear.  Cylindri¬ 
cal  lenses  could  be  used  to  focus  the  light  preferentially  in  one  dimen¬ 
sion  to  provide  a  strip  of  illumination  at  much  lower  total  power.  How¬ 
ever,  cylindrical  lenses  are  difficult  to  grind  and  are  often  just  slabs 
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of  glass  rod.  Some  applications  would  warrant  the  design  of  special 
lenses  and  spatial  filtering,  but  for  most  installations  the  power  loss  is 
not  prohibitive  and  is  compensated  for  by  ease  of  design  and  availability 
of  components. 


3.5  PARTI CLE-ON-A-DISK  STUDY 

Another  method  for  determining  depth  of  field  effects  was  to  utilize 
a  glass  bead  fixed  to  a  glass  disk.  The  disk  was  attached  to  a  variable  - 
speed  electric  motor,  and  an  accurate  number  of  revolutions  could  be 
obtained  by  using  a  strobe  unit  and  timer. 

Unfortunately,  some  potential  problems  exist  with  this  method 
which  cannot  be  totally  eliminated.  First  is  the  selection  of  a  glass 
disk.  Since  small  particulates  and  scratches  will  show  up  as  lower- 
channel  noise,  they  must  be  minimized.  The  particle  was  attached  to  the 
glass  by  a  weak  solution  of  sugar  and  water.  Second,  distributions  ob¬ 
tained  in  this  manner  are  not  statistically  representative  of  even  a  class 
of  particles  centered  about  the  size  of  interest.  This  is  caused  by  repe¬ 
tition  of  the  same  particle,  which  will  track  through  approximately  the 
same  portion  of  the  array  each  revolution. 

Such  a  study  was  performed  using  a  3  80 -Mm  glass  bead.  Figure  22 
is  a  composite  distribution  obtained  by  summing  each  channel  count 
over  the  entire  depth  of  field.  This  is,  in  effect,  a  simulation  of  a  uni¬ 
form  distribution.  From  0  to  1/4  farfields  the  distribution  was  centered 
on  450  Mm  with  some  counts  in  the  first  channel.  The  low -channel  counts 
were  attributable  to  dust  and  surface  imperfections  and  will  not  be  in¬ 
cluded  in  the  composite  distribution.  From  1/4  farfield  on  out  to  1-1/4 
farfields  the  average  measured  size  increased,  finally  stabilizing  at  an’ 
indicated  size  of  625  Mm.  This  correlates  with  the  first  diffraction  ring 
beyond  the  particle  edge.  Data  were  taken  using  a  5-mw  He-Ne  laser 
with  a  0  of  about  0.55.  While  this  intensity  ratio  is  not  sufficient  for 
the  first  diffraction  ring  to  trigger  by  itself,  it  does  drop  the  intensity 
to  a  point  where  added  noise  may  cause  false  triggering.  This  conclu¬ 
sion  can  be  further  justified  by  the  gradual  increase  in  measured  size 
with  distance.  It  is  this  latter  effect  which  is  perhaps  the  most  serious 
limitation  of  the  FOS  as  it  is  presently  designed. 
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SECTION  IV 
APPLICATIONS 


4.1  ROCKET  NOZZLE  TEST 

Tests  were  performed  on  a  small  rocket  injector  with  water  simu¬ 
lating  both  the  fuel  and  the  oxidizer.  Normally  NTO/MMH  or  NTO/UDMH 
in  the  ratio  of  about  2:1  is  used.  The  injector  was  similar  to  those  used 
in  the  Apollo  Command  Service  Module  (CSM)  and  Lunar  Module  (LM) 
reaction  control  thrusters  with  a  design  flow  rate  of  approximately  0.35 
lbm/sec.  The  hole  pattern  is  shown  in  Fig.  23.  Figure  24  illustrates 
the  detailed  dimensions  of  a  simple  hole  doublet. 

This  particular  injector  was  chosen  since  data  existed  concerning 
the  particle  distribution  as  measured  using  an  in-line  holocamera 
(Ref.  6  and  Fig.  25).  Water  was  injected  at  75  psig  into  both  the  oxi¬ 
dizer  and  fuel  manifolds.  This  resulted  in  an  average  discharge  veloc¬ 
ity  of  2.  7  x  10^  cm/sec  at  the  nozzle  face. 

Data  were  taken  6.  1  cm  downstream  of  the  rocket  face  and  2.  54  cm 
off  the  axis.  Masking  tubes  provided  mechanical  aperturing  of  the  spray. 
An  upstream  slit  width  of  7.  49  cm  was  maintained  while  the  downstream 
side  was  cut  away  to  eliminate  entrained  water  from  settling  on  the 
optics.  Figure  26  shows  the  masking  tubes.  Splitter  rings  were  turned 
on  the  tubes  to  help  prevent  droplets  from  flowing  along  the  surface  and 
off  the  ends  into  the  sampling  area.  Some  breakup  of  droplets  caused 
by  the  edge  was  expected.  Since  the  slit  width  was  approximately  10 
times  the  depth  of  field  of  the  smallest  particle,  the  effect  of  breakup  on 
the  result  should  be  small.  The  end  effects  were  therefore  assumed  to 
be  negligible. 

The  bin  size  was  selected  as  25  ptm,  and  a  30 -percent  threshold 
was  used.  The  illumination  source  was  a  5-mw  He-Ne  laser  with  a 
beam  diameter  at  the  exit  of  about  0.  65  mm  at  the  1/e2  point.  Beam 
divergence  was  1.7  milliradians.  This  provided  an  illumination  spot 
diameter  at  the  sample  area  of  approximately  1.  50  mm.  This  is 
approximately  twice  the  array  image  or  maximum  measurable  particle 
size. 

Table  II  summarizes  the  results  of  the  FOS  data  taken  at  the  same 
conditions  as  the  holography  test.  An  explanation  of  the  symbols  and 
format  follows  the  table.  For  comparison  a  constant  velocity  for  all 
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particles  was  assumed  which  was  applied  to  the  measured  particle  flux 
to  obtain  the  densities.  An  average  density  of  24.  56/ cm^  was  measured 
over  a  period  of  7.  44  sec  as  compared  with  33  ±  5.  2/cm^  at  one  instant 
in  time  for  the  holography  (as  based  on  2  7  particles).  The  particle  flux 
as  measured  by  the  FOS  is  shown  in  Fig,  27.  The  log  of  the  flux  is 
plotted,  as  there  usually  is  a  larger  number  of  particles  in  the  small 
channels.  However,  these  particles  contribute  little  to  the  overall 
volumetric  average  diameter.  Also  shown  are  the  geometric  sizes  of 
the  injector  holes  (see  Fig.  24). 

A  Nukiyama-Tanasawa  plot  of  both  the  FOS  and  the  holography  data 
was  made  in  which  the  two  are  compared  (see  Fig.  28  and  Ref.  7).  The 
Nukiyama-Tanasawa  expression  derived  for  jet  impingement,  following 
Ingebo,  is  given  by 


log  F  =  log 


where 


Volume  contribution  of  size  D  particles 
Total  volume  of  all  particles 

=  bin  size,  and 

=  volume  mean  diameter. 


In  general,  Dmp  is  given  by 


D  -M' 

mp  l?ni  DipJ 


A  plot  of  log  F  versus  D  has  a  slope  given  by 
Dgo  =  —1.7 /Slope 


In  Fig.  28  a  volumetric  mean  diameter  was  calculated  directly  from 
both  FOS  and  holographic  data.  This  was  used  to  generate  the  lines 
passing  through  the  data. 

The  data  were  corrected  for  two  reasons.  First,  although  the  same 
bin  size  was  used,  the  method  of  sorting  the  particles  by  the  FOS  was 
basically  different  from  the  manual  reduction  method  used  with  the  holo¬ 
graphic  analysis.  An  approximate  1/2 -bin  shift  was  introduced  which 
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which  was  adjusted  at  the  lower  end.  Second,  the  probability  for  count¬ 
ing  the  large  particles  is  very  low.  In  fact,  if  a  total  of  400  particles 
is  recorded,  as  was  the  case  for  the  hologram,  the  probability  of  not 
counting  a  particle  in  the  higher  channels  becomes  greater  than  the 
probability  of  obtaining  a  count.  Therefore,  the  data  obtained  by  the 
FOS  was  again  modified  to  exclude  these  points  for  comparison  pur¬ 
poses  only. 

Agreement  is  in  general  good;  however,  severe  deviations  are 
noted  at  both  the  high  and  low  ends  of  the  range.  As  this  effect  is 
common  to  both  sets  of  data,  it  will  appear  that  the  irregularity  is  a 
function  of  the  fluid -particle  dynamics  rather  than  a  system  error. 

The  time  average  D30  =  293.4  ^m  of  the  FOS  was  based  on  16,061 
points.  The  spatial  average  D30  =  270  pm  of  the  holographic  system 
was  based  on  400  points.  The  difference  is  approximately  8.  3  percent 
or  0.  936  b. 

Plume  profiles  of  various  quantities  can  also  be  readily  obtained. 
With  the  aperture  tubes  in  place  a  scan  was  made  over  the  radius  of  the 
plume  at  a  distance  of  23.  5  cm  downstream.  The  aperture  width  was 
fixed  at  1.  19  cm,  thereby  sampling  an  area  of  approximately  2.  5  mm2. 
Figure  29  shows  plots  of  diameter  as  a  function  of  plume  radius.  Note 
the  decrease  in  size  as  the  edge  of  the  plume  is  reached.  Hologram 
shadowgraphs  near  the  face  indicate  this  trend.  Figure  30  plots  parti¬ 
cle  flux  and  density  as  a  function  of  radius.  This,  too,  clearly  falls 
off  in  a  monotonic  fashion.  Particle  flux  measurements  are  exact  in¬ 
sofar  as  this  is  the  measured  quantity.  However,  the  inferred  density 
should  be  taken  with  some  caution,  as  the  discharge  velocity  at  the 
face  was  used  in  its  calculation.  A  decrease  in  particle  velocity  with 
increasing  radius  and  downstream  distance  would  be  expected.  This 
velocity,  in  general,  is  size  dependent.  An  independent  velocity  mea¬ 
surement  could  be  used  to  produce  densities  if  desired. 


4.2  DUST  EROSION  TUNNEL 

The  FOS  was  installed  in  the  PWT  Dust  Erosion  Tunnel  to  investi¬ 
gate  system  performance  at  high  velocities  of  reasonably  well-known 
particles. 

The  flow  is  obtained  from  a  continuous  operation  arc-heated  wind 
tunnel  equipped  to  conduct  erosion  tests  of  materials  under  high-velocity 
and  high-enthalpy  conditions  (see  Fig.  31).  The  tunnel,  which  has  a 
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maximum  run  time  of  five  minutes,  is  capable  of  using  particles  from 
50  to  1000  pm  with  velocities  up  to  6500  ft/sec  depending  on  size.  It 
utilizes  high-pressure  air  heated  in  a  5-MW  arc  heater. 

Particles  are  injected  upstream  of  the  nozzle  throat  of  a  low- 
expansion-rate,  hypersonic  nozzle,  where  they  are  aerodynamically 
accelerated.  The  test  section  consists  of  a  test  cabin  and  a  multiple- 
mount  model  injection  system  (see  Fig.  32).  An  exhaust  connection  is 
provided  through  a  diffuser  to  the  Engine  Test  Facility  (ETF)  exhaust 
plant.  The  installation  is  shown  in  Fig.  33. 


4.3  MICROSCOPE  COMPARISON 

A  comparison  of  grits  prepared  for  use  in  the  Dust  Erosion  Tunnel 
(DET)  was  made.  Nominal  6  50 -Pm  grit  was  sized  by  the  FOS  utilizing 
the  aspirator  and  dispenser.  A  sample  was  also  sized  by  conventional 
microscope  techniques.  Figure  34  shows  a  direct  comparison  of  the 
two  distributions.  Both  distributions  show  a  very  large  number  of 
particles  in  the  5-  to  100-pm  range  in  addition  to  those  greater  than 
650  jtim.  The  FOS  did  not  include  particles  below  16  pm,  which  may 
account  for  the  lower  number  of  particles  registered  in  this  bin. 

Figure  35  shows  the  same  distributions  normalized  to  the  total  number 
of  particles  in  each  distribution. 

A  notable  departure  is  found  in  the  200-  to  650-pm  region.  This 
is  a  direct  result  of  the  use  of  magnesium  oxide  particles  and  the 
method  used  in  characterizing  the  diameter.  Figure  36  is  a  photo¬ 
micrograph  of  the  particles.  Note  the  irregular  nature  of  the  crystaline 
particles  and  the  presence  of  small  fragmentation  particles.  The  FOS 
sizing  is  based  on  a  projected  length.  Random  orientation  of  the  parti¬ 
cles,  as  verified  by  double-pulsed  holography,  produces  the  wide  vari¬ 
ation  in  measured  sizes.  The  microscope  count  was  based  on  the 
maximum  linear  dimension.  The  average  ratio  of  maximum -to¬ 
rn  inimum  linear  dimension  is  1.4.  Furthermore,  a  close  inspection 
under  the  microscope  indicates  that,  due  to  their  geometry,  the  partic¬ 
ulates  orient  themselves  preferentially  on  the  surfaces.  In  general, 
they  seek  a  stable  orientation,  thereby  introducing  a  bias  in  the  distri¬ 
bution  toward  the  larger  indicated  size. 
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SECTION  V 
APPARATUS 


5.1  THEORY  OF  OPERATION 

The  fiber  optics  system  has  already  been  described  in  some  detail. 

It  consists  of  an  imaging  system  of  sufficient  numerical  aperture  to  pass 
at  least  the  first  minimum  of  the  smallest  particle  pattern.  This  portion 
must  be  specifically  designed  for  particular  application.  Next  is  the 
fiber  array,  extending  for  0.  844  cm  and  composed  of  thirty-two  264-jum 
fiber  optics.  The  electronic  detection  system  consists  of  32  PM  modules, 
an  input  sensor  strobe  module,  and  a  size  encoder  (see  Fig.  37).  Ex¬ 
ternal  power  supplies  provide  -6,  +12,  and  +5  volts  and  an  adjustable  high- 
voltage  supply  for  the  photomultipliers .  An  additional  clock  and  memory 
inhibit  were  added  to  provide  timed  input  for  rate  determination. 

5.2  PHOTOMULTIPLIER  DETECTION  MODULE 

As  shown  in  Fig.  37,  light  is  channeled  from  the  linear  array  to  the 
corresponding  931A  PM  tubes  located  on  printed  circuit  boards  in  a 
modular  chassis  rack.  Interlocking,  light  baffled  tubes  allow  the  boards 
to  be  removed  for  servicing  and  still  allow  a  close  interface  between 
fiber  ends  and  the  PM  tubes.  The  current  from  the  tubes,  nominally 
100  ma,  is  connected  to  the  input  of  a  702  current -to -voltage  converter 
producing  a  signal  level  of  about  3  volts.  The  output  of  the  702  is  sampled 
and  is  inverted  by  a  741  operational  amplifier  to  provide  a  reference  volt¬ 
age.  A  low -pass  input  filter  allows  power  supply  ripple  and  1/F  noise  to 
pass  to  later  null  with  similar  noise  in  the  signal. 

The  low-pass  3-db  point  can  be  increased  up  to  almost  0.  5  MHz, 
where  the  741  begins  to  roll  off.  The  signal  and  reference  outputs  are 
summed  at  the  input  of  a  710  comparator.  The  comparator  changes 
state  if  the  signal  drops  below  the  predetermined  threshold.  Some 
feedback  is  provided  for  clean  crossover.  This  is  supplied  by  voltage 
division  of  the  positive  output  and  typically  is  50  mv. 

The  comparator  output  is  inverted  by  a  7402  NOR  inverter,  the  out¬ 
put  of  which  is  through  a  pair  of  diodes  (IN914).  One  diode  from  each 
of  the  32  modules  leads  to  a  common  output  to  form  an  expandable  NAND 
function  to  sense  the  particle  transit.  The  other  diode  output  is  inverted 
and  sets  a  latch.  The  Q-latch  output  connects  to  the  size  determination 
encoder,  while  the  Q  output  drives  a  lamp  through  a  2N3643  transistor 
for  static  tests  of  the  individual  boards.  The  latches  are  reset  by  pulses 
from  the  input  sensor  and  strobe  generator  logic  module. 
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5.3  INPUT  SENSOR  AND  STROBE  GENERATOR 

The  common  outputs  from  the  PM  board  modules  form  the  NAND 
Gate  input  (Fig.  37),  which  is  at  a  logic  '0'  during  particle  shadow 
transit.  A  DT932  quad -input  OR  gate  provides  a  logic  'l1  for  the  transit 
time.  This  output  is  inverted,  and  the  trailing  edge  of  the  pulse  sets  an 
S-R  flip-flop  (7474).  The  leading  edge  of  the  7474  output  triggers  a 
9601  one-shot,  providing  a  delay  to  allow  for  crossover  noise.  The  out¬ 
put  triggers  a  second  one-shot  to  provide  a  strobe  pulse.  The  last  one- 
shot  provides  a  reset  pulse  for  the  PM  module  latches  after  a  suitable 
delay.  The  7474  output  also  is  used  to  provide  an  inhibit  pulse  back  to 
the  PM  modules  to  prevent  a  second  particle  from  entering  until  the  first 
has  been  processed. 

5.4  SIZE  DETERMINATION  ENCODER 

The  size  encoder  (Fig.  37)  provides  a  parallel  binary  representation 
of  the  memory  address  to  the  particle  data  handling  system.  The  out¬ 
puts  from  the  PM  modules  are  summed  in  five  adding  levels  by  a  combi¬ 
nation  of  full  adders.  The  basic  circuitry  is  standard.  However,  in 
addition  to  the  30  size  channels  possible  are  the  conditions  of  single  and 
reject  and  oversize  or  double-end  reject.  These  channels  are  0  and  31, 
respectively,  and  correspond  to  triggering  of  one  or  both  of  the  end 
modules  during  a  particle  transit. 

5.5  PARTICLE  DATA  SYSTEM 

Up  to  this  point,  the  electronic  processing  has  been  described  in 
some  detail.  The  data  handling  system,  however,  contributes  little  to 
the  information-gathering  process  and  merely  provides  data  storage, 
display,  and  handling  capabilities  and  will  be  only  briefly  discussed 
(see  Fig.  38). 

The  data  handling  system  provides  storage  for  32  channels  of  infor¬ 
mation,  30  of  which  correspond  to  size  data;  the  remaining  two  are 
used  for  end-reject  data  as  previously  described.  Each  of  these 
channels  may  be  read  out  by  the  Nixie  tube  display  or  on  the  cathode 
ray  tube  (CRT).  They  are  sequentially  updated  at  the  frame  rate. 

(One  full  frame  is  equivalent  to  a  complete  reading  of  the  32  memory 
units. )  The  memory  may  be  reset  manually  or  automatically  by  selec¬ 
tion  of  a  frame  reset  rate.  Memory  transfer  is  controlled  manually, 
but  may  be  controlled  by  an  external  source.  This  simply  controls  the 
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acceptance  of  data  into  the  memory  units.  In  conjunction  with  the 
memory  transfer  switch  an  external  scaler-timer  is  enabled  to  give 
elapsed  counting  time. 

The  CRT  display  can  be  scaled  by  up  to  a  factor  of  100,  calibrated, 
and  positioned  from  the  panel.  An  analog  signal  is  also  provided, 
which  is  a  direct  result  of  the  digital-to- analog  conversion  necessary 
for  the  vertical  sweep. 

BCD  output  is  also  provided,  together  with  a  print  command  pulse, 
to  drive  a  printer.  The  channel  number  and  the  counts  in  that  channel 
are  printed  out  sequentially  at  a  rate  which  can  be  determined  through 
hardwiring  to  the  timing  oscillator. 


5.6  REJECT  PROCESSOR  AND  GATED  STROBE 

Information  from  the  particle -sizing  system  includes  five  binary  lines 
for  particle  sizes,  two  lines  for  end  fibers,  and  a  strobe  line  (Fig.  37). 
The  logic  for  acceptance  or  rejection  of  the  size  determination  is  per¬ 
formed  within  the  particle  data  system.  It  is  at  this  point  that  a  valid 
strobe  pulse  is  generated  to  allow  the  size  determination  to  address  a 
memory  unit.  Assuming  that  no  end  fibers  have  been  occulted,  the  BCD 
input  will  be  used  to  address  the  appropriate  memory  and  the  strobe 
pulse  will  index  the  count  contained  by  one.  An  invalid  or  reject  strobe, 
if  generated  by  one  or  both  fibers  being  occulted,  will  cause  the  appro¬ 
priate  address,  0  or  31,  and  index  it  by  one. 

Additional  gating  is  provided  to  drive  counters  to  record  separately 
right-  or  left-end  fiber  rejects.  This  information  might  be  useful  for 
highly  nonuniform  flows. 


SECTION  VI 

CONCLUDING  REMARKS 


The  FOS  offers  a  versatile  means  of  particle  sizing  in  flow -field 
situations.  It  has  real-time  capability  and  is  compatible  with  most 
digital  data  handling  systems.  Data  reduction  programs  are  straight¬ 
forward  and  proven.  Through  the  use  of  curve  fits  and  prerun  calibra¬ 
tion  a  minimal  amount  of  computing  capability  is  required  for  an  on-line 
monitoring  system. 


27 


AEDC-TR-73-1 1 1 


The  major  drawback  of  the  system  is  that  caused  by  poor  electronic 
depth-of -field  cutoff.  This  effect  can  be  reduced  by  means  of  mechanical 
aperturing  as  was  done  in  the  water  spray  tests.  Other  approaches 
exist  which  could  reduce  the  effects  of  noise  coupled  with  the  out-of- 
focus  particle  transits.  A  crossed  beam  coincidence  would  reduce 
these  spurious  effects. 

The  crossed  beam  approach  would  rely  on  the  coincidence  of  scat¬ 
tered  light  and  particle  shadow  transit.  The  volume  would  be  defined 
by  that  volume  about  the  object  plane  of  the  imaging  system  in  common 
with  an  off-axis  collection  system  and  an  apertured  photomultiplier  tube. 
Such  a  system  would  allow  localized  measurements  to  be  made  within 
flow  fields  whose  extent  along  the  optical  axis  is  on  the  order  of  or  less 
than  the  cutoff  distance. 
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Fig.  4  Normalized  Mean  Size  versus  Channel  for  Various  Threshold  Values 
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Fig.  5  Normalized  Array  Width  versus  Channel  Number  for  Various  Thresholds 
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Fig,  6  In-Focus  Particle  on  Array 
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Fig.  8  Particle  on  Linear  Fiber  Array 
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Fig.  13  Intensity  Distribution  of  a  20-/nm  Particle 
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Fig.  15  Traverse  System 
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Fig.  16  Particle  Dispenser 


AEDC-TR  -73-1  1 1 


Cutoff  Length,  cm 


AEDC-TR-73-1 1 1 


Particle  Diameter,  p.m 

Fig.  17  Cutoff  Length  versus  Particle  Diameter 
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Fig.  18  Particle  in  Gaussian  Beam 
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Fig.  21  Focused  Gaussian  Beam 
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Fig.  25  Rocket  Injector  Spray  Test 
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Fig.  27  Particle  Flux  versus  Size 
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Fig.  29  Mean  Diameter  versus  Radius 
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Fig.  30  Particle  Flux  and  Density  versus  Radius 
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Fig.  31  Layout  of  Dust  Erosion  Tunnel 
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Fig.  32  Optical  Arrangement  through  Test  Cabin 
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Fig.  33  Fiber  Optics  Particle  Sizer  in  the  Dust  Erosion  Tunnel 
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Fig.  38  Particle  Data  System 
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TABLE  I 

DATA  REDUCTION  PROGRAM 


FORTRAN  IV  G  LEVEL  20  MAIN  DATE  -  73131  11/22/13 


C 

SET00216  OC  TOOO  1/13/72  B85252  HISTOGRAM  PLOT  1  REVISED) 

0031 

DIMENSION  TITLE!  201  *0(301  «W(30)  , PI 30 >  ,0N  131)  , ELOO) ,  A <30 > ,F ! 30) , 
1  GI30I >DEN( 301 fR( 30) ,GF(30) 

0002 

PEAL  *8  ENO«EN<  30) «EN31 ,SEN 

0003 

DIMENSION  8UFF  12178) 

0004 

CALL  CLASS  (BUFF, 8712) 

0005 

WRITE ( 6, 2007) 

0006 

10 

REA0I5, 10C0, END-900) TITLE 

0007 

REAOI 5, 1002)8, T.C.E.V.ELX, THRESH, START 

0008 

JSTR- START 

0009 

REAOI 5, 10C2)ENO,EN,EN31 

0010 

CALL  PDF!A,THRESH,D,W,P) 

0011 

DM! 11*01 1 )-• 5*8 

0012 

DO  20  J-2,30 

0013 

20 

0M(J)-.5*!D(J )+D( J-l) ) 

0014 

ONI  311*0! 30)4,5*8 

0015 

00  30  J-l, 30 

COlb 

EL! J) -AMIN  11 .5*C*0( J) **2/E,ELX) 

0017 

AIJI-ELI J)*W( J) 

0018 

FIJ)-EN(J)/IPI J)*EL(J)*T) 

C019 

G1 J)-FI J )*( DM( J+1)-DM( J)) 

0020 

OENI J)-GIJ)/V 

0021 

30 

CONTINUE 

0022 

S-0. 

0023 

00  40  J-l, 30 

0024 

IF(J.LT.JSTRIGD  TO  39 

0025 

S-S»GIJI*DI Jl**3 

0026 

39 

CONTINUE 

0027 

RIJI-S 

0028 

40 

GF(JI-G(J)/ia*OI J)**2I 

0029 

DO  52  J-l, 30 

0030 

RIJ)*R(J)/S 

0031 

IFIGFI J).GT.O.)GO  TO  51 

0032 

GF1J)— 1.E50 

0033 

GO  TO  52 

0034 

51 

GF! J)=AL0G10IGF! J)/S) 

0035 

52 

CONTINUE 

0036 

SOEN-O. 

0037 

SG=0. 

0038 

SEN-0. DO 

0039 

010-0. 

0040 

C20-0. 

0041 

0  30-0. 

0042  ' 

DO  60  J-JSTR ,30 

0043 

SOEN=SDEN+OfcNI J 1 

0044 

SG=SG+G(J) 

0045 

SEN-SEN+EN 1 J ) 

0046 

010-010+G! J)*0(  Jl 

0047 

D20-D20+G! J  )*D( J )**2 

0048 

60 

D30-D30+G! J>*0!  J)**3 

0049 

D10-010/SG 

0050 

020- SORT ( 020/SG ) 

C051 

D30-I 030/SG )**• 3333333 

0052 

WRITE (6, 2000) TITLE 

0053 

WRITE  16,2001) 

0054 

NO-ENO 
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TABLE  I  (Continued) 


FORTRAN  IV  G  LEVEL  20  MAIN  OATE  »  73131  11/22/13 


0055 

N31-EN31 

0056 

WRITE! 6, 3002IN0 ,N31 ,B  ,T*C ,E , V ,ELX*THRE$H, JSTR 

0057 

WRITEI6.2003) 

0058 

DO  70  J=l,30 

0059 

N*ENI J) 

0060 

70 

WRITEI6,2002IJ,N,D< J) ,G ( J I ,DENl J) ,EL I J) ,W( Jl , A< Jl t R( J) . GFU I > FI J ) 

0061 

N-SEN 

C062 

WRITE! 6,2005  IN. SG, SOEN 

0063 

WRITE (6, 2006) DIO, 020,030 

0064 

WRITEI6, 20071 

0065 

CALL  GRABIOM, G, TITLE, •  D  * , *  J  *) 

0066 

CALL  GRAFID.F, TITLE,'  0  F  *1 

0067 

GO  TO  10 

CQ68 

900 

WRITE  16, 2004 ) 

0069 

CALL  KLASS 

0070 

STOP 

0071 

1000 

FORMAT! 20A4I 

0072 

1002 

FORMAT (6E 12* 01 

0073 

2000 

FORMAT! 1H1.19X.20A4) 

0074 

2001 

FORMAT!  IH0.T4.6HS INGLE, T1  3, 6H DOUBLE, T 109, 8HST ART I NG/ 

1  T4.6HREJECT.T13, 6HR EJECT ,722, 
21HB, T34, 1HT ,T46 , 1HC ,T  58,6HLAMB0A ,T70 , 1HV ,T82, 4HLMAX ,T94,6HTHRESH, 
3T109, 7HCHANNEL) 

0075 

2002 

FORMAT! 21 9, 1P9E 12*41 

0076 

3002 

FORMAT 1219, IP 7E 12*4,1 121 

0077 

2003 

FORMAT! 1H0,T9,1HI ,T17 , 1HN , T22 ,1 HD.T34 ,1H J,T46 .7HDENS ITV ,T58 , 1HL, 

IT 70, 1HW.T82, 1HA , T94 , 1HR.T 106 ,6H LOG! F ) ,T1 17.9HDIST  FUNC) 

0078 

2004 

FORMAT! 'OTHE  END') 

0079 

2005 

FDRMATI 9H0  TOTALS, 19 , 12X , IP2E12.4I 

0060 

2006 

FORMAT! 1H0.T 14, 5HD10  *, 1PE12.4/T14.5H020  -, E12.4/T14, 5HD30  -,E12.4 
1) 

FORMAT! 1HAI 

0081 

20C7 

0082 

END 
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TABLE  I  (Continued) 


FORTPAN  IV  G  LEVEL  20  PDF  DATE  -  73131 


0001 

SUBROUTINE  POFI B,C ,D,W,E> 

0002 

DIMENSION  D(30)iM(30)«E(30)fN(30) 

0003 

T=l.-C 

0004 

DO  1  J*=l,30 

0005 

DIJ>*0. 

0006 

W(J)=0. 

0007 

1 

NIJI-0 

0008 

DO  3  L=1 » 161 

0009 

X«.1*IL-.5I 

0010 

DO  2  M-1,321 

0011 

DH=.05*(H-.5I 

0012 

AP-DH+X 

0013 

JP«AP 

0014 

EP=AP-JP 

0015 

CALL  LUNEYIDH.EP) 

0016 

IFIEP.GT.T)JP«JP*1 

0017 

IFI JP.GT.15IG0  TO  2 

0018 

AM-DH-X 

0019 

JM=AM 

0020 

I  F(  AM  .LT  .  0.  )  JH=  JN—  1 

0021 

FH=AM-JM 

0022 

CALL  LUNEYIDH.EM) 

0023 

IFIEN.CT.T1JN-JN+1 

0024 

IFI JM.GT. 151G0  TO  2 

0025 

JFIJM.LT.-15IG0  TO  2 

0026 

J-JP+JH 

0027 

IFIJ.LT.1IG0  TO  2 

0028 

NIJ)«NI  JIH 

0029 

Dl J )-0l J  )*0H 

0030 

2 

CONTINUE 

0031 

3 

CONTINUE 

0032 

DO  10  J=l,30 

0033 

10 

DIJ)=DIJ)/N(J) 

0034 

DO  30  L~ 1 > 161 

0935 

X=.1*!L-.5I 

0036 

DO  20  K-1.30 

0037 

OH“DI  X 1 

0038 

AP-DH+X 

0039 

JP-AP 

0040 

EP=AP-JP 

0041 

CALL  LUNE YIDH«EP) 

0042 

IFIEP.GT.TIJP-JP41 

0043 

I F(  JP.GT.  15 )GO  TO  20 

0044 

AH-DH-X 

0045 

JM*AM 

0046 

IFIAM.LT. 0.) JH=  JM-1 

0047 

EM-AH-JM 

0048 

CALL  LUNEVIOHfEMI 

0049 

IFIEN.GT.TI JH-JMU 

0050 

I  FI JM.GT . 15)G0  TO  20 

C051 

IFIJM.LT.-15IG0  TO  20 

0052 

J* JP* JM 

0053 

I  FI J.NE.K1G0  TO  20 

0054 

WIJI-WI J»*l. 

0055 

20 

CONTINUE 

0056 

30 

CONTINUE 
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TABLE  I  (Continued) 


FOP TP AN  IV  G  LEVEL  20 


POF 


OATE  «  T3L31 


0057 

DO  4  J»1 i 30 

0058 

D(JI=2.*B*0I J) 

0059 

J) 

0060 

A 

FIJI* .02*B*B*NI  J 1 

0061 

PFTURN 

0062 

ENO 

FORTRAN  IV  G  LEVEL  20  LUNEV  DATE  -  T3131 


0001 

SUBROUTINE  LUNEYIDH.EPI 

0002 

REAL  PI/3.141593/ 

0003 

INTEGER  K/0/ 

0004 

D«2.*DH 

0005 

IFIEP.LT. DJGO  TO  10 

0006 

EP=D*D 

0007 

RETURN 

0008 

10 

CONTINUE 

0009 

A»2.*EP 

0010 

A-A/IA-l.-DI 

0011 

CT1=«1.-A*IEP-D> 

0012 

CT2«1.-A*(EP- 1.1/0 

0013 

STI»S0RT(1.-CT1*CT1J 

0014 

ST2-SQRTU.-CT2*CT2J 

0015 

T1«ATAN?( STltCTlI 

0016 

T2-ATAN2IST2.CT2) 

0017 

EP*( T1-ST1*CT1»D*0*I T2-ST2*CT2I 1 /PI 

0018 

RFTURN 

0019 

ENO 
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TABLE  I  (Continued) 

FORTRAN  IV  G  LEVEL  20  GRAB  OATE  »  73131  11/22/13 


0001 

SUBROUTINE  GRAB (A, B, TITLE «LBX,LBY) 

0002 

DIMENSION  A<  31 J  tB<30)  (TITLE 1 20)  ,XI120)  (VI 120)  (RI2)(SI2) 

0003 

R 1 1)“ AMINII 0. ( AI1) ) 

0004 

P I2)>AI31 ) 

0005 

SID-1.E50 

0006 

SI2)-0. 

0007 

DO  10  J=l«  30 

0008 

K-441 J-l) 

0009 

XIK+1  l=AI'J) 

0010 

XIK«2)-AIJ) 

0011 

XlKOI-AIJ+l) 

0012 

XIK*4I*AI J*l) 

0013 

Y 1 K* 1 ) =0. 

0014 

VIM-2I-BI J) 

0015 

YIK+31-BI Jl 

0016 

VIK*4)~0. 

0017 

IFIBI J1.LE.0.IG0  TO  1C 

0018 

SI1I-AHIN1ISI1) (BIJ)) 

0019 

SI2)*AMAX1ISI2),B1J)I 

0020 

10 

CONTINUE 

0021 

CALL  CALCMP (l.,l»,0,2) 

0022 

CALL  SCAL EI2(l(Ril,7(LXd X, JX.MX.NX) 

0023 

CALL  SCALEI2(1(S(0,C(LY(IV,JY,MY(NY) 

0024 

MY=1 

0025 

NY«1 

0026 

J Y=I V46 

0027 

CALL  AXISIl.,1. ,0.,-13. ,LX,IX,JX,NX,NX,4,LBX) 

0028 

CALL  AX) SI  1.(1. (90. ,8.0,LY, IYf JY,MY,NY,4,LBY) 

0029 

CALL  HASHI1.,13.(LX,IX( JX,KX,NX(1..8. (LY(IY, JY,MY(NY> 

0030 

CALL  SYMBOL  11. (9. 25, .21, TITLE (0..72) 

0031 

CALL  LINOI 0(0(0,1,X,1. ,13. ,LX(I X, JX(HX, NX,XN,XX ,0X( 

1  1 , Y, 1 . ,  8.,LY.IY,JY, MY,NY,YN,YX,OY) 

0032 

YMN=10.**VN 

0033 

DO  20  J  =  l, 120 

0034 

20 

Y 1 J 1 = ANAX 1 1 Y( J 1 , YMN) 

0035 

CALL  LINEI 120, 0,0,1, X, LX, l.Y.LY) 

0036 

RETURN 

0037 

ENO 
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TABLE  I  (Concluded) 


FORTRAN  IV  G  LEVEL  20 


GRAF  DATE  -  73131 


0001 

SUBROUTINE  GRAFIA.B.TITLE.LBX.LBY) 

0002 

DIMENSION  A!  301  .81301  ,TI TLEI 20)  .XI 30 >  ,Y  1 30)  »R (21 .  SI2> 

0003 

RIll-AMINlIO.fAIll) 

0004 

R!2)>AI 30) 

0005 

S(1)=1.E50 

0006 

SC  2)»0. 

0007 

DO  10  J-1,30 

0008 

X  C  J)= AC  J 1 

0009 

YC  J)3BC J) 

0010 

IFCBC JI.LE.O.JGO  TO  10 

0011 

SI1)>AMIN1(S!1I,BM|) 

0012 

SC2)=AMAX1( SC2) .BCJI ) 

0013 

10 

CONTINUE 

0014 

CALL  CALCHPI1..1..0.2) 

0015 

CALL  SCALE ( 2.1.R.1.7.LX.IX.JX.MX.NX) 

0016 

CALL  SCALE! 2. l.S. 0.0, LY.I Y.JY.MY.NY) 

0017 

MY-1 

0018 

N  Y=1 

0019 

J Y=IY*6 

0020 

CALL  AXISIl.,1. .0..-13. .LX. IX. JX.MX, NX.4.LBX) 

0021 

CALL  AXIS! 1..1. ,90. .8.0.LY.I Y, JY.MY.NY.4.LBY) 

0022 

CALL  HASHI1«,13..LX.IX.JX,MX.NX,1.,B>  »LY ,IY,JY,HY,NYI 

0023 

CALL  SYMBOL! 1., 9. 25,. 21, TITLE, 0., 72) 

0024 

CALL  LINO!0,0.0.1.X.1..13.,LX,IX.JX.MX,NX,XN,XX,DX. 

1  I.Y.I.,  8..LY,IY.JY.NY,NY,YN.YX.0Y) 

0025 

YMN=10.**YN 

0026 

DO  20  J-1.30 

0027 

20 

YC  J)«AMAX1!  YC J),YMN) 

0028 

CALL  LINEI30. O.O.l.X, LX. l.Y.LV) 

0029 

RETURN 

0030 

END 

11/22/13 
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TABLE  II 

SAMPLE  COMPUTER  PRINTOUT 


AEDC -TR-73-111 


SYMBOLS 


B 
T 
C 

Lambda 
V 
L 

max 
Single  reject 
Double  reject 
I 

D 
N 
J 

Density 
L 
W 
A 
R 

Log  (F) 

Totals  are  self-explanatory. 

Dio  =  Simple  average 

D20  =  Area  weighted  average 

DgQ  =  Volume  weighted  average 

"Thres,  "  if  printed,  is  the  electronic  threshold  value  and 
has  a  bearing  on  the  exact  Dj  used  in  making  calculations. 

Log  (F)  when  plotted  versus  diameter  yields  a  Nukiyama- 
Tanasawa  plot. 

The  graphs  are  plots  of  particle  flux.  The  first  is  based  on  predust 
conditions.  A  total  particle  flux  of  less  than  20  particles/ cm2-sec  was 
noted.  The  particles  were  about  50  Mm  in  diameter.  Note:  particles 
below  approximately  one  bin  size  are  not  counted,  which  may  lead  to 
erroneous  conclusions.  This  can  be  a  very  important  consideration. 


=  Bin  size  (e.  g. ,  3157  jum) 

=  Run  time  in  seconds 
=  Constant  for  depth  of  field 
=  Wavelength  of  source 

=  Assumed  particle  velocity  for  density  calculation 
=  Maximum  allowable  depth  of  field 

=  Number  of  times  one  end  fiber  was  shadowed 
=  Number  of  times  both  end  fibers  were  shadowed 
=  Bin  number 
=  Particle  diameter 
=  Number  of  particles  measured  in  bin 
=  Particle  flux  (number/  cm^ -sec-bin) 

=  Particle  density  (number  / cm^  -bin) 

=  Length  of  optical  path 
=  Effective  array  width 
=  Effective  array  area 
=  Cumulative  volume  fraction 

=  ar 

AD  D5 
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TABLE  II  (Concluded) 


#4 

X-0.V-1.1«Z*2.4  9/21/71 

ROCKET  INJECTOR  HOLOGRAPHY  COMPARISON 

B  = 

2.50006-03 

T  * 

7.4445E 

oc 

c  * 

2.5000E  00 

LAMBDA  ■ 

6.32806-05 

V  ' 

2.6960E 

03 

LHAX  = 

7.4930E-01 

SINGLE  REJECT  - 

9610. 

D0U8LE  REJECT  = 

612. 

1 

N 

D 

J 

RHO 

L 

W 

A 

R 

LOGIFI 

1 

2549. 

2. 5000E-03 

3.6492E 

04 

1.3536E  01 

1.2346E-01 

7.6000E-02 

9, 3829E-03 

2.  1939E-04 

1.1954E 

01 

2 

2038. 

5. 0000E-03 

7. 5422 E 

03 

2.7976E  00 

4.9384E-01 

7.3500E-02 

3.6297E-02 

5.8215E-04 

1.0667E 

01 

3 

1752. 

7. 5000E-03 

4.4237E 

03 

1.6408E  00 

7.4930F-01 

7 . 1000E-02 

5. 3200E-02 

1.3002E-03 

1.0083E 

01 

4 

1768. 

1.0000E-02 

4.6270E 

03 

1.7163E  00 

7.4930E-01 

6.8500E-02 

5 . 1327E-02 

3.0806E-03 

9.8526E 

00 

5 

1203. 

1.2500E-02 

3.2676E 

03 

1.2120E  00 

7.4930E-01 

6.6000E-02 

4.9454E— 02 

5.5362E-03 

9.5077E 

00 

6 

888. 

1. 5000E-02 

2.5070E 

03 

9.2988E-01 

7. 4930  E-01 

6.3500E-02 

4. 7580E-O2 

8.7918E-03 

9.2342E 

00 

7 

821. 

1.7500E-02 

2.4128E 

03 

8.9496E-01 

7.4930E-01 

6.1000E-02 

4.5707E-02 

1.3767E-02 

9.0837E 

00 

8 

805. 

2. 000OE-02 

2.4665E 

03 

5.  1502E-01 

7.4930E-01 

5.8500E-02 

4.3834E-02 

2. 1361E-02 

8.9774E 

00 

9 

718. 

2*  2500E-02 

2.2985E 

03 

8.5256E-01 

7.4930E-01 

5.6000E-02 

4.1961E-02 

3. 14356-02 

8.8443E 

00 

10 

613. 

2.5000E-02 

2.0541E 

03 

7.61 90E—01 

7. 4930 E-01 

5.3500E— 02 

4.0088 E-02 

4.3784E-02 

8.7040E 

00 

11 

460. 

2.7500E-02 

1. 6  WOE 

03 

5.9976E-01 

7.4930E-01 

5.1000E-02 

3.8214E-02 

5.67  23E-02 

8.5173E 

00 

12 

397. 

3.0000E-02 

1.4674E 

03 

5.4430E-01 

7 • 4930E-01 

4.8500E-02 

3.6341E— 02 

7.1968E-02 

8.3996E 

00 

13 

251. 

3. 2500E— 02 

9.7819E 

02 

3.6283E— 01 

7. 49306—01 

4.6000E-02 

3.4468E-02 

8.4888E— 02 

8.1539E 

00 

14 

190. 

3. 5000E-02 

7.8302E 

02 

2.9044E-01 

7.4930E-01 

4.3500E-02 

3.2595E-02 

9.7806E-02 

7.  992 9E 

00 

15 

168. 

3.7500E— 02 

7.3457E 

02 

2.7247E-01 

7.4930E-01 

4.1000E-02 

3.0721 E— 02 

1 . 1 27  IE— 01 

7.9052E 

00 

16 

350. 

4.0000E-02 

1.6297E 

03 

6.0450E-01 

7.49*06-01 

3.8500E-02 

2.8848E-02 

1.5284E— 01 

8.1953E 

00 

17 

112. 

4.2500E-02 

5.5773E 

02 

2.0687E-01 

7.4930E-01 

3.6000E-02 

2.6975F— 02 

1.6932E-01 

7.6769E 

00 

18 

141. 

4. 5000E-02 

7.5454E 

02 

2.7987E-01 

7.4930E-01 

3.3500E-02 

2.5 102E-02 

1.9577E-01 

7.7585E 

00 

19 

82. 

4.7500E-02 

4.7420E 

02 

1.7589E-01 

7.4930E— 01 

3 • 1000E-02 

2. 322BE— 02 

2. 1533E-01 

7.509SE 

00 

20 

84. 

5.0000E-02 

5.2838E 

02 

1.9599E— 01 

7.4930E-01 

2.8500E-02 

2.1355E-02 

2.4074E— 01 

7.5123E 

00 

21 

56. 

5.2500E-02 

3.8612E 

02 

1.4322E-01 

7. 4930 E-01 

2.6000E-02 

1.9482E-02 

2.6224E-01 

7.3337E 

00 

22 

36. 

5.5000E-02 

2.7463E 

02 

1.0186E-01 

7.4930E— 01 

2 .3  500E-02 

1.7609E-02 

2 .7982E-dl 

7.1453E 

00 

23 

20. 

5.7500E— 02 

1.7C73E 

02 

6. 3329E-02 

7.4930E-01 

2 . 1000E-02 

1.5735E-02 

2.9231E-01 

6.9003E 

00 

24 

117. 

6. OOOOE— 02 

1.1333E 

03 

4.2C5*E-01 

7.4930E-01 

1 . 3  500E-C2 

1 . 3862E-02 

3. 9654E-01 

7.6855E 

00 

25 

259. 

6.2500E-02 

2.9C19E 

03 

1.0764E  00 

7.4930E-01 

1.6000E-02 

1.1989E— 02 

6.5914E-01 

8.0582E 

00 

26 

87. 

6.5000E-02 

1.1553E 

03 

4.2852E-01 

7.4930E-01 

1 . 3500E-02 

1.01 16E-02 

7.8122E-01 

7.6241E 

00 

27 

66  • 

6. 7500E-02 

1.0756E 

Oi 

3. 9897E-01 

7.4930E-01 

1 . 1000E-02 

8.2423E— 03 

9.0851E— 01 

7.5603E 

00 

28 

25. 

7.0000E-02 

5.2727E 

0? 

1.9557E-01 

7.49 30 E-01 

0. 5000 E -03 

6. 3690E-03 

9.7809E-01 

7.2191F 

00 

29 

5. 

7. 2500E-02 

1.4939E 

02 

5.5413E-02 

7.4930E-01 

6.0000E-03 

4.4958E-03 

1. OOOOE  00 

6.6409E 

00 

30 

0. 

7. 5000E-02 

0.0 

0.0 

7. 49 30 E-01 

3.5000E-03 

2.6225E-03 

1. OOOOE  00 

-1. OOOOE 

50 

TOTAL  N 

030 

TOT  KHO 

16061. 

3.1222E-02 

3.1673E  01 
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TABLE  til 

INTENSITY  DISTRIBUTION  PROGRAM 


FORTRAN  IV  G  LEVEL  20 


HAIN  OATE  =  73131'  11/23/54 


0001 

0002 

0003 

00C4 

0005 

0006 

0007 

3308 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0330 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0942 

0043 

0044 

0045 

0046 

0047 

0348 

0049 

0050 

0051 

0052 

0053 

0054 


C  SET00244  OC  TOOD  06/13/72  BB5252  HT  BENTLEY 

C  INTENSITY  OISTP IBUTION  3RD  REVISION  DIRECT  METHOD 

REAL* 8  RAI200I.ARE1200) »BRE( 2001 1 AAGI200 1 • BAG 1 200) «TENI200 1 

CALL  CLASS! 'SET00244* I 

P 1=3. 141593 

T0L-1.E-4 

S«44. 

DZ-O. 

EM=1. 

EL-5. 08 
F-22. 

0- .002 

W-4.765E-5 

Z-S+OZ 

2P-EM*S 

CAYH-PI/W 

CAY-2. *CAYH 

ELH- . 5*EL 

A-CAYH*! l./ZP-l./PI 
B-A  *C AYH/Z 

IF(4BS(B*ELH*ELH).LE..01IB-0. 

NP-15 
NT-20000 
WRITE! 9,10051 

WRITE ( 9. 10061 S >02 >EM|EL  >F .0. NR 
WRITE  I 9i 1007 ) 

WRITE! 9» 10061 WiZtZP ,CAY»A,B»NT 

NP1-NR-1 

NT1-NT-1 

0H-.5*0 

BC— CA YH*0/Z 

RT-1. 5*EN*0 

DP-RT/NR1 

0T-ELH/NT1 

T80-C AY*RT*ELH/ZP 

TB1-BC*FLH 

TAT-A*ELH*ELH. 

T8T-B*ELH*ELH 

WRITE (6> 10061 TBOtTBl • TAT. TBT 
DR- 1.0-4 
DO  1  J-l.NR 
PA! JI-I J-1»*0R 
ARE! J)=0.00 
BPEI J)-0.00 
AAGIJ 1-0.00 
1  BAG!  J )  -0.00 

P A< 151-30.0-4 
Cl-CAY/ZP 
DO  5  K  =  1  *  NT 
T-l K-l )*DT 
TS-T*T 
AR-A*TS 
SNA-S IN ( ARI 
CSA-COS! AR1 
AR=B*TS 
SNB-S INIARI 
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AEDC-TR-73-11 1 


FORTRAN 

0055 
0056 
COST 
0058 
0059 
C060 
0061 
0062 
0063 
0064 
0065 
0066 
0067 
0068 
0069 
OOTO 
0071 
0072 
0073 
0074 
0075 
0  076 
0077 
0078 
0079 
0080 
0081 
C082 
0083 
0084 
0085 
0086 
0087 
0088 
0089 
0090 
0091 
0092 
0093 
0094 
0095 
CO  96 
0097 
0098 
0099 

0100 

0101 


TABLE  III  (Continued) 


IV  G  LEVEL  20  MAIN  DATE  -  73131  11/23/S4 

CSB-COS(Aft) 

N*1 

X«8C*T 

CALL  BESJ!X,N«8J,T0L,  JER) 

IF! JER.EQ.OIGO  TO  2 

WRITE 1 6,90031 X»BJ»TOL ,N, JER 

2  BJ1-BJ 
N«0 

C2“C1*T 
DO  4  J»1,NR 
X»C2*RA( J| 

CALL  RESJI X,N,B J,TOL,  JER  I 
IFI JER.EQ.OIGO  TO  3 
HRITE(6,9000>X,BJ,T0L,N,JER 

3  BJO-BJ 
BBJ»8J0*B41 
TBJ-T*BJO 

ARE( JI-AREl JI*TBJ«CSA 
BRE! J>  «BREI J1 ♦8BJ*SNB 
AAGIJfAAGI  JI+TBJ*SNA 
BAG! JI'BAGI Jt*B8J*CSB 

4  CONTINUE 

5  CONTINUE 
TAN=C1*C1 
C3»0H*0T 

DO  6  Jal.NR 
ARE! JI»DT*ARECJJ 
RPE ( J I »C3*8RE ( J I 
A  AG! J1»0T*AAGIJ1 
BAG! JlaC3*BAGIJI 

TENI Jt>TAN*l I ARE! JI-BRE ( J) )**2+lAAG!JI»BAG(Jl 1**2) 

6  CONTINUE 

CALL  SCALEX1NR.RA,'  R2  »l  PLOT 

CALL  GRAFIRA, ARE, 'REAL  A',6) 

CALL  GRAFIRA, BRE, 'REAL  B',61 
CALL  GRAFIRA.AAG,' IMAGINARY  A', lit 
CALL  GRAFIRA, BAG, 'IMAGINARY  8', lit 

CALL  GRAFIRA, TEN, 'INTENSITY', 9J  PLOT 

WRITE  I  6, 10041 

CALL  CLASS! 'SET00244'! 

STOP 

1004  FORMAT! •  THE  END' I 

1005  FORMAT! 4X, 1HS ,1 IX ■ 2H0Z , 10X.1HM. 11X,1HL,11X,1HF,11X, 1HD, 13X, 2HNR I 

1006  FORMAT! 1P6E12. 4,1 8t 

1007  FORMAT!  4X,6HL  AMBDA,6X  ,1HZ  ,  11'X,2HZP,10X,1HK,11X,  1HA  ,1 IX,  1H8, 13X.2HN 
IT  I 

9000  FORMAT! 1P3E 14.6,21141 
END 
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TABLE  III  (Continued) 


FORTRAN  IV  G  LEVEL  20 


GRAF 


OATE  -  73131  11/23/54 


0001 

0072 

0003 

0004 

0005 

0006 

0007 

0006 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0016 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 


SUBROUTINE  GR  AF  (  X,  Y,  LABY  ,  NL) 

DIMENSION  X|1),YI1)  .LABYI  1>  ,LA6  (201 
INTEGER  K/2/ 

CALL  SCALE  (N,K,Y,5,4,LY, IT. JY.MY.NY) 

CALL  CALCMP  II. ,1. ,0.2) 

RFWIND  9 
P  EAOI 9, 210 )LAB 

CALL  SYMBOL  I  1.,  10. 25,.  19,  LAB,  0.  ,80) 

PFAOI  9, 2101LAB 

CALL  SYMBOL! 1. ,10.00, .19, LAB, 0. , 80) 

REAOI 9.210ILAB 

CALL  SYMBOL! 1., 9.50,. 19,LAB,0., 80) 

R E AO! 9, 2 10 ) LAB 

CALL  SYMBOL!  l.,9.  25,.  19,  LAB,  0.  ,80) 

CALL  HASH! 1., 13., LX, I X, JX, MX, NX, 1., 8., LY, IV, JY.MY.NY) 
CALL  AXIS!1.,1.,0. ,-13.,LX,IX,JX,MX,NX,4,LABX) 

CALL  AXIS! l.,l.,90.,8.,LY,IY, JY ,MY, NY ,NL  .LABY ) 

CALL  L I NOIN ,0,0 ,K,X,1.,13. ,LX,IX, JX.HX, NX, 01 ,02 ,03, 

1  K.Y.I.,  8.,LY,IY,JY,HY, NY ,04 , OV ,06) 

RETURN 

c 

ENTRY  SC1LEXIN1 , X.LBX ) 

N=N1 

LABX'LBX 

CALL  SCALE! N,K, X,5, 7.LX.I X, JX,MX,NX) 

RETURN 

210  FORMAT! 20A4) 

ENO 


FORTRAN  IV  G  IFVEL  20  BESJ  DATE  -  73131  11/23/54 


0001 

SUBROUTINE  BESJ! X,N,B J,0, IER) 

BESJ 

042 

0002 

PEAL  PI/3.141593/ 

0003 

8J-.0 

BESJ 

044 

0004 

IEP'0 

0005 

IFIN ) 10, 20,20 

BESJ 

045 

0006 

10 

IER-1 

BESJ 

046 

0007 

RETURN 

BESJ 

047 

0008 

20 

IF! X)30, 25, 31 

0009 

25 

IFIN. EU.0)BJS1. 

0010 

RETURN 

0011 

30 

IEP-2 

0012 

PFTURN 

BESJ 

050 

0013 

31 

I F<  X.GT . 25. )GQ  TO  310 

0014 

IF(X-15. 132,32, 34 

0015 

32 

NTCST=20. *10. *X— X**  2/3 

BESJ 

052 

0016 

GO  TO  36 

BESJ 

053 

0017 

34 

NTEST=93.*X/2. 

BESJ 

054 

0018 

36 

IFIN-NTESTI 40,38, 38 

BESJ 

055 

0019 

38 

I  ER*4 

BESJ 

056 

0020 

RETURN 

BESJ 

057 
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TABLE  III  (Concluded) 


FOB  "PAN 

IV  G  LEVEL 

20  BESJ  DATE  -  73131 

11/23/54 

0021 

40 

IFP  =  0 

BESJ 

05  B 

0022 

N1=N+1 

BESJ 

059 

0523 

PPPEV=.0 

BESJ 

060 

0024 

IFIX-5. >50, 60,60 

BESJ 

064 

0025 

50 

MA*X+6. 

BESJ 

065 

5526 

GO  TO  70 

BESJ 

066 

0027 

60 

HA*1.4*X*60./X 

BESJ 

067 

0028 

70 

MB*N+lFlX(X>/4+2 

BESJ 

06B 

0029 

MIERO*MAXOMA,MB> 

BESJ 

069 

0030 

mmax=ntest 

BESJ 

073 

0031 

100 

no  190  M*M2ER0,P<NAX,3 

BESJ 

074 

0532 

FM1=1.0E-2B 

BESJ 

07B 

0333 

FM=.0 

BESJ 

079 

0034 

ALPHA*. 0 

BESJ 

OBO 

0035 

IF(M-(M/2)*2I120,11C,120 

BESJ 

081 

0036 

110 

JT*-1 

BESJ 

0B2 

0037 

GO  TO  130 

BESJ 

0B3 

0038 

120 

JT-1 

BESJ 

0B4 

0039 

130 

M2=M-2 

BESJ 

0B5 

0040 

CO  160  K-1.M2 

BESJ 

0B6 

0041 

MK=M-K 

BESJ 

087 

0042 

PMK=2.»FL0AT<MKI*FH1/X-FM 

BESJ 

088 

0043 

FH*FM1 

BESJ 

0B9 

0044 

FM1*BMK 

BESJ 

090 

0045 

I F( MK— N- 11150, 140, 150 

BESJ 

091 

0046 

140 

B J*BHK 

BESJ 

092 

0047 

150 

JT*-JT 

BESJ 

093 

0048 

S  =  1«-JT 

BESJ 

094 

0049 

160 

alpha*alpha+bmk*s 

BESJ 

095 

0050 

BPK*2  .*FH  1/X-FM 

BESJ 

096 

0051 

1F(N ) 183, 173,180 

BESJ 

097 

0052 

170 

BJ=BMK 

BESJ 

098 

0053 

180 

ALPH*=ALPHA*B1K 

BESJ 

099 

0054 

BJ=BJ/ALbHA 

BESJ 

100 

0555 

I F(ABS(BJ-BPREV>-ABS(D*BJ> >200,200, 190 

BESJ 

101 

0056 

190 

BPPEV=8J 

BESJ 

102 

0057 

IFP  =  3 

BESJ 

103 

0058 

2  00 

RETURN 

BESJ 

104 

0359 

310 

BJ=SQPTI2./(PI*X>>*C0S<X-.5*PI*(N-.5>  > 

0060 

PFTUPN 

0061 

END 

BESJ 

10  5 
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graphic  technique  to  determine  the  dimensions  and  numbers  of  particles  moving  in  a 
flow  field.  The  system  is  digital  in  nature.  Particles  pass  through  a  collimated 
laser  beam  and  are  imaged  onto  a  linear  array  by  a  coaxial  lens.  The  array  is  com¬ 
posed  of  the  exposed  ends  of  a  fiber  optics  bundle  which  serves  as  a  "link"  between 
the  array  plane  and  the  sensing  photo-detector  modules.  Being  an  imaging  device, 
it  can  measure  a  wide  range  of  particle  sizes  through  the  proper  selection  of  optics. 
Sizes  ranging  from  2  to  1500  p  m  have  been  measured  in  the  course  of  this  project. 
Comparisons  of  holographic  data  of  a  liquid  rocket  injector  and  of  water  spray  nozzles 
are  made  with  the  fiber  optics  system  . 
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